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Motivated by atom optics experiments, we investigate a new class of fidelity functions 
describing the reconstruction of quantum states by time-reversal operations as Mos.{t) = 



me 



iH2t/2 iHit/2 -iH2t/2 ~iHit/2 



IV') I • We show that the decay of Mua. is quartic in time at short 



times, and that it freezes well above the ergodic value at long times, when H2 — Hi is not too large. 
The long-time saturation value of Mua contains easily extractable information on the strength of 
decoherence in these systems. 



PACS numbers: 05.45.Mt, 03.65.Yz 



I. INTRODUCTION 

When subjected to external noisy fields, quantum me- 
chanical wavefunctions lose memory of their phase. As 
a fundamentally important consequence of this decoher- 
ence process, pairs of partially scattered waves no longer 
interfere, and the dynamics follows the Liouville time- 
evolution of classical densities A somehow similar 
situation occurs when one time-evolves an initial super- 
position = Catpa of many eigenmodes of the 
Hamiltonian Hi governing the time evolution, with in- 
commensurate eigenfrequencies £„. In this case, for each 
pair of components (a,/?), the relative phase (e^ — ep)t 
becomes pseudo-random, which washes out partial wave 
interferences. This dephasing process, however, differs 
from decoherence in a fundamental way that it can in 
principle be undone by an appropriate time-inversion. As 
a matter of fact, echo experiments are able to reverse the 
sign of the Hamiltonian, Hi —Hi, by means of effec- 
tive changes of coordinate axes induced by electromag- 
netic pulses 0. When this operation is performed after 
an evolution time t, one expects the initial wavefunc- 
tion to be reconstructed at 2t, regardless of its spread 
over eigenmodes. Imperfections in the pulse sequence or 
unavoidable couplings to external uncontrolled degrees 
of freedom result instead in an imperfect time-inversion. 
Hi — > —Hj = —Hi — Y,, and therefore the Loschmidt 
echo dMl i^l (we set ^ = 1) 



ruLit) = (V-l 



with 



,iH2t-iHit 



(la) 
(lb) 



gives a better description of the fidelity with which the 
experiment reconstructs the initial state. Echo experi- 
ments in nuclear magnetic resonance 0, QiCfuantum op- 
tics atomic [E43|, condensed matter Q, microwave 
cavities Q, and elastodynamics [l3| have demonstrated 
that Ml(<) remains sizable for times significantly longer 
than the dephasing time. The decay of Mi^{t) allows one 
to extract information on irreversible decoherence pro- 
cesses induced by S. 



In experiments with cold atoms the Loschmidt echo, 
Ml, can be extracted from interference fringes of Ramsey 
spectroscopy [ll|. There, an effectively two-level atom is 
initially prepared in a state |1) ® ['(/'), where |1) and |2) 
denote the two internal atomic states, and stands for 
the spatial component of the initial state. First, the atom 
is irradiated with a microwave frequency field with energy 
chosen to change the atomic state into an equiprobable 
superposition of |1) ® |'0) and |2) ® \'4>). Such a field 
is referred to as a 7r/2 pulse. The atom is then let to 
evolve in an optical trap for a time t during which the 
l)-component of the state evolves under a spatial Hamil- 
tonian Hi, while the |2}-component under H2. After that 
another 7r/2 pulse is applied to the atom and the proba- 
bility P2 for the atom to be found in the internal state |2) 
is measured. It turns out that this probability is essen- 
tially determined by the Loschmidt echo amplitude, ttil. 
In practice however one works not with a pure initial 
state IV'), but with a thermal mixture of initial states. 
The echo amplitude ttil from each of these states con- 
tributes to P2 with a different, effectively random phase, 
which in turn reduces the fringe contrast in a Ramsey ex- 
periment. As a result, the 7r/2-7r/2 pulse sequence proves 
inefficient in measuring the Loschmidt echo for large en- 
sembles of thermally populated states. 

In order to overcome this difficulty Davidson and col- 
laborators implemented a novelpulse sequence in their 
echo spectroscopy experiments [a, 01 

MDa(i) = |mDa(i)P , with (2a) 
mDa(t) = {ilAe'"-"^e'"'"^e-'"^"^e-'"'"^\ilj) . (2b) 

The corresponding pulse sequence consisted of three 
short pulses, 7r/2-7r-7r/2, separated by two time intervals 
of equal duration t/2, after which P2 was measured. The 
TT pulse swaps the population of the internal states |1) 
and 1 2). The probability P2 is then determined by the 
amplitude rriDa, and each individual state of the ther- 
mal ensemble contributes to P2 with the same phase. 
Thus, the 7r/2-7r-7r/2 pulse sequence allows one to mea- 
sure the echo in Eq. ([2]) even for ensembles of more than 
10^ of thermally populated states, as in the experiments 



2 



of Refs. @,0|. 

It is clear from the definitions given by Eqs. ([T]) and ^ 
that mathematically Moa is not the same quantity as the 
Loschmidt echo Ml . Even though some significant differ- 
ences between Moa and Mj^ have been previously envis- 
aged in the literature they have never been systematically 
studied. It is the purpose of this article to fill in this gap 
by comparing the two quantities both analytically and 
numerically. Below, we show that Moa differs from the 
Loschmidt echo A/l in the two important respects that (i) 
its short-time decay is quartic and not quadratic in time, 
and (ii) for not too strong perturbation S = H2 — Hi, 
AfDa saturates at a perturbation-dependent value, well 
above the ergodic saturation of Ml (00) ^ N^^ at the 
inverse Hilbert space size. Fidelity freezes have been re- 
ported for Loschmidt echoes with off-diagonal perturba- 
tions with a zero time average , phase-space displace- 
ment perturbations [l3| . and more recently for initially 
pure states coupled to complex environments [T^ . how- 
ever the freeze we report here has a different physical 
origin. We note in particular that it persists for i — > 00. 
The long-time saturation of Moa allows to extract the 
strength of the fields in E more easily than by fitting 
decay curves of conventional echoes, over not precisely 
defined time intervals. Moreover, the absence of decay 
arising from S in this pulse sequence makes it straightfor- 
ward to extract decoherence rates because, assuming that 
the pulse sequence is perfect, any decay in experimentally 
obtained data for MDa(i) would come exclusively from 
the coupling of the system to external degrees of freedom, 
not included in our theory. Given the superb experimen- 
tal control that modern echo experiments have on their 
pulse sequence, this novel echo spectroscopy has therefore 
the potential to deliver precious, previously unattainable 
information on the dominant sources of decoherence in 
trapped cold atomic gases. 

II. SHORT-TIME DECAY 

There has been a large number of analytical and nu- 
merical investigations of the Loschmidt echo and some of 
its offsprings (l5l.[l6|. Most, if not all approaches assume 
a small perturbation, i.e. |S| <C 1^^1,2 1 for an appropri- 
ate operator norm. As but one consequence, the largest 
energy scale is the energy bandwidth B, which to lead- 
ing order is the same for Hi and H2- For short times, 
t <Si B~^, Mi^{t) is easily calculated by expanding the 
propagators in Eqs. ([T]), and keeping the leading order 
contributions. One obtains 

ML(i) ~ 1 - (aLi)' , (3) 

where 

^2 = Sl = Fi-F2. (4) 

Thus, the short-time decay of the Loschmidt echo is 
quadratic (itI . [18| . with a rate given by the dispersion 



(Tl of the perturbation operator El evaluated over the 
initial state. 

The same procedure can be applied to Moa, where it 
however gives 

MDa(0 ^ 1 - (^Dai)' , (5) 

with the decay rate (TDa given by 

= {^A^lM - (^ISDal^)' , SDa = |[i^l,i^2] ■ 

(6) 

Two things are remarkable here. First, the short-time de- 
cay of Afoa is quartic in t, and thus slower than the decay 
of A/l. Second, its rate is determined by the commutator 
of the unperturbed and perturbed Hamiltonians. 

III. LONG-TIME SATURATION 

The analysis of the long-time behavior of A/l and Mbs, 
starts by diagonalizing the unperturbed and perturbed 
Hamiltonian operators, Hi = "^^^ Eu\u)(u\ and H2 = 
X]-u and expanding the initial state in the basis 

of the unperturbed Hamiltonian, = X^u^ul'") [3- 
The resulting expression for the echo is then averaged 
over time to yield the mean saturation value. In the case 
of the Loschmidt echo the time-averaged saturation is 
given by 

Ml,oo= <c^"\cu'\'{u\v){v\u'){u'\v){v\u"). (7) 

The next step is to average this expression over a ran- 
dom ensemble of coefficient c„ for the initial state, such 
that c*Cu' = N~^du.u'- (Hereinafter, an overline de- 
notes the averaging over an ensemble of random initial 
states.) Here N is the effective size of the Hilbert space 
(the number of eigenstates of Hi 2 comprising the initial 
state). To leading order in 1/N, one uses c*c„"|ck'P = 
c*Cm" • \cu' |2 = N^'^Su,u" to obtain the ergodic saturation 
value 

M^^j^. (8) 

Using the same procedure, one can calculate the long- 
time saturation value of AfDa- At the level of the echo 
amplitude mDa, one gets 

"iDa,oo = ^ c*^Cu'{u\v){v\u'){u'\v){v\u') (9a) 

u.u' ,v 

= 4ElH«)r- (9b) 



One then uses an approximation | (ulu) | ~ | | with 
"i Ts" 

= p{Eu — Ey) a function of only the energy differ- 
ence between the two states. Replacing one of the sums 
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in Eq. ^ by an integral over the energy difference be- 
tween the two states scaled by the mean level spacing 
A = B/N, we can write 



dE 
A 



pHe) 



(10) 



This expression relates the long-time saturation of 
to the energy spreading of eigenfunctions of Hi over those 
of H2 as measured by p{E). It is known for a large variety 
of quantum chaotic systems that, in the regime A <C 
r/A ^ B, this spreading has a Lorentzian shape 



piEu - E„) 



A 



r/2 



TT {Eu - E^Y + (r/2)2 



fii) 



with a spreading width given by the golden rule, T ~ 
(Tl/A [ill, see Eq. Q for the definition of o-£. We thus 
obtain 



rUDa 



A_ 



fl2) 



Eqs. (fTTj) and ^ predict an average saturation value 
A^Da,cx) above the ergodic saturation for N < (B/ttT)'^. 
The width T of the Lorentzian ([TT|) increases with |S|, 
and the ergodic saturation, Eq.®, is recovered when T > 
B/ttTVI/^ thus 



AV - 

ttE j ' N 



(13) 



We note that the Lorentzian spreading of Eq. ([TT|) is re- 
placed by more complicated, system-dependent spiked 
structures in dynamical systems with mixed or regular 
dynamics, for which it is accordingly impossible to draw 
general conclusions. We stress, however, that Eq. ((TU]) 
remains valid even in that case. 

Eq. is the main result of this paper. This new 
long-time fidelity saturation originates from the specific 
sequence of time-evolutions in Moa, giving the long-time 
behavior of the latter as an energy integral over the 
squared average overlap of eigenstates \u) of Hi 

over the eigenstates \v) of H2- For completeness we next 
comment on the intermediate regime, between the short- 
time quartic decay and the long-time saturation. 



IV. INTERMEDIATE ASYMPTOTIC DECAY 

We briefly sketch a semiclassical analysis of Afoa in the 
intermediate regime between the short-time decay and 
the long-time saturation. We follow the lines of Ref. (20| 
to show that AfDa and A/l have the same behavior in 
that regime. 

In the semiclassical approximation the time evolution 
of under Hj, j = 1,2, is given by 



(re 



7(r'^r,t) 



e*^^-^(r'IV') . (14) 




Figure 1: (Color online) Trajectories of the unperturbed (71) 
and perturbed (72, 72, and 72) systems together with the 
associated action differences {AS, AS", and AS"). 



Here, the sum goes over all classical paths 7 connecting 
r' and r in time t, Sj,-y = Sj^^{r,r' ,t) is the action along 
7, Dj,^ = (2ni)-'^^^\dct{d^Sj,^/drdr')\^/^e-''"'^--'^^ with 
Morse index Vj^j counting the number of conjugatepoints 
on 7, and d is the dimensionality of the system |2l| . The 
semiclassical Loschmidt echo amplitude is obtained by 
inserting Eq. in Eq. (jlb[) . The resulting expres- 

sion contains three spatial integrals over r, r', and r" 
with a double sum over trajectories 71 (r' — !■ r,t) and 
72 (r" r,t) corresponding to the Hamiltonians Hi 
and H2 respectively. The standard analysis of this ex- 
pression involves three steps (i) One assumes that 
{r\ip) is localized about a point Tq and evaluates the inte- 
grals over r' and r" by stationary phase approximations. 
This reduces the set of paths 71 and 72 to those start- 
ing at ro, see Fig. [T] (ii) Noting that the double sum 
over trajectories contains rapidly oscillating phase fac- 



tors exp[i(S'i_^j — 52.72)] J so that only pairs of correlated 
paths 7i and 72 contribute to mDa, one employs the di- 
agonal approximation (72 ^^_7i) to reduce ttil to a sum 
over a single path 71 . Ref. |22| , building up on ideas first 
expressed in Ref. [1^, justified this step by the shadowing 
theorem, (iii) Finally, one uses the fact that iDi,^! P is 
the Jacobian of a transformation between final positions 
r and initial momenta p on paths 71 . This allows one to 
change the integration variable from r to p to get 



mL(t) = (27r) 



dp e'^^ I (p 



(15) 



Here AS* = AS'(ro,p,t) = Si^^-^ — S'2,72 the difference 
between the action of an unperturbed trajectory 71 leav- 
ing the point Tq with a momentum p and traveling for 
time t and the action of the corresponding perturbed tra- 
jectory 72 ~ 7i . Following the same procedure one finds 



nT-Da(t) = (27r)" 



dp e 



i(A5'-AS") 



|(p|V')| 



(16) 



where AS" = AS'(ro, p, i/2) and A5"' AS'(f,p,i/2) 
with (f , p) being the phase space point on 71 at time 
t/2, see Fig.[Tl In other words, AS" (AS"') is the action 



4 




Figure 2: (Color online) Average echo decay for the kicked 
rotator model with Ki = 57, N = 8192, and K2 - Ki = 
5 • 10'^ 1.2 ■ 10"*, 2.1 • 10"* and 5 ■ 10"* (Ml, red solid 
lines from top to bottom), and K2 — Ki = 5 • 10"^, 9 • 10"''', 
1.2 • 10"*, 1.6 • 10"*, 2.1 ■ 10"*, 3.1 • 10"*, and 5 ■ 10"* (Afoa, 
black dashed lines, from top to bottom). Curves are averages 
over 500 initial states. 

difference between the first (second) half of the unper- 
turbed trajectory 71 and the corresponding perturbed 
trajectory 72 (72)- This is sketched in Fig. [TJ 

Once averaged over an ensemble of initial states, both 
Mi^ and Afoa satisfy 



ML,Da(i) ^ mL,Da(t) (17) 

+(2^)-2'^ [dp f dp'i(piv;)n(p'iv^)i', 

J J Up 

where the integral over p' is restricted to a volume Op 
around p, such that two trajectories starting from the 
same spatial point with momenta p and p' G ftp stay 
"close" in phase space during time t. The first term in 
the right-hand side of Eq. ((T7|) is evaluated using the cen- 
tral limit theorem, exp(zAS') ~ exp(— AS'2/2) ~ e~'"*/^ 
and exp[i{AS' - AS")] ~ exp[-(AS^ + AS^)/2] ~ 
g-r(t/2+t/2)/2 ^ e-r*/^ where T is defined in Eq. 1^ 
as the width of the local density of states. For A/Da, 
we neglect correlations between AS" and AS"', which is 
justified by the fast decay of correlations along chaotic 
classical trajectories. The second term in Eq. P7)) is de- 
termined by the measure of the set Op and in chaotic 
systems decays as e""^* with A being the average Lya- 
punov exponent of the underlying classical system |2C| . 
Therefore, the intermediate time decay of Afoa is the 
same as that of A/l [Sflj, [23|, i-e. 

MlM MdI(I) - e-* """[^'^1 . (18) 

This exponential time decay continues until the echo 
reaches the saturation plateau given by Eq. (|13p . 

Ref. reported some saturation of Af^a for ultra- 
cold atoms inside optical traps. However, at this stage, a 



direct comparison of these experiments with our theory 
does not seem feasible, because they explore completely 
different time regimes. Indeed, the echo spectroscopy 
experiments of Refs. [1,01 are concerned with short times 
corresponding to no more than 3-4 oscillations/bounces 
of an atom in the trap. In contrast, the semiclassical 
derivation of the exponential decay, Eq. (|18p . and the 
RMT analysis of the fidelity freeze, Eq. p^ . are only 
valid for times much longer than the average free flight 
time. 



V. NUMERICAL STUDY 

We confirm our analytical results with some numerical 
data. Our simulations are based on the kicked rotator 
model with dimensionless Hamiltonian 

i?l,2 = y +/^l,2C0s£^(5(t-71T). (19) 

n 

For large enough kicking strength, /vi.2T > 7, the dy- 
namics is fully chaotic with a Lyapunov exponent A = 
ln[i4ri,2T/2]. We quantize this Hamiltonian on a torus, 
and accordingly consider discrete values pi = 2ttI /N and 
xi = 2ttI/N, I = 1, ...A^, giving an effective Planck's 
constant h^s = l/N. Both echoes Mi^{n) and AfDa('T-) 
are computed for discrete times t ~ nr, with the 
kicking period r, using the unitary Floquet operators 
Ui^2 = eyip[—ip^/2hes]exTp[—iKi^2COSx/hc{f] for single- 
kick time-evolutions. The bandwidth is _B = 27r and ac- 
cordingly A = 2'k/N . The eigenstates of U2 spread over 
those of Ui according to Eq. ((TTI) with F oc {6K N)'^, 
with SK = K2 - Ki [Hi- Together with Eq. we 
thus expect a long-time saturation of A^Da at a value 

M^^{dK^NY\ (20) 

for SK^N^ < 1. 

Fig. [5] shows the time decay of the echoes, Mi^{t) 
shown as red curves and A'/Da(i) as black curves, aver- 
aged over an ensemble of randomly chosen initial states. 
For equal values of the perturbation strength, both A/l 
and A^Da display an exponential time decay governed 
by the same decay rate, providing a clear support for 
Eq. . The Loschmidt echo decay saturates at a value 
^ A^"^ in agreement with Eq. (|5]). The freeze of A/Da 
occurs at a value that decreases with increasing per- 
turbation strength until it reaches ergodic saturation at 
N~^. We confirm in Fig. [3] that the numerically ob- 
served perturbation-dependent saturation of A^Da follows 
Eq. (|20p . Once plotted as a function of SKN^/^, satura- 
tion data for N e [256,8192] and SK G [A ■ 10"^ 0.052] 
nicely fall on top of one another until they deviate be- 
cause they have different ergodic saturation, N~^. More- 
over, in the regime of validity A ^ 6K <C S of Eq. pT|) . 
one has AfDa.oo oc (SKN^^^)'' with an exponent b ~ 3.8 
close to the prediction b — A from Eq. . We note that 



5 




Figure 3: (Color online) Long-time saturation value of Afoa 
for A'l — 57 and A'^ = 256 (black up triangles), 512 (red cir- 
cles), 1024 (green squares), 2048 (blue down triangles), 4096 
(violet diamonds) and 8192 (cyan left triangles). Main panel; 
rescaled data confirming the analytical prediction of Eq. (|20|) . 
The straight black line indicates a slope of oc Inset: 
raw data as function of the difference of dimensionless kicking 
strengths 5K = K2 — Ki, with -R'1,2 = -^'1,2/1". 



b is larger for data with larger Hilbert space size iV, where 
the fitting range is larger - and the fit is accordingly more 
accurate - because saturation occurs at larger values of 
SKN^/^. We also checked numerically that the initial 
decay of A/Da is quartic and not quadratic in time. Our 
numerical simulations thus fully confirm the theoretical 



predictions derived above. 



VI. CONCLUSIONS 

Our analysis of the new fidelity function Mua @, 01 
shows that it significantly differs from the Loschmidt 
echo in two important respects: (i) the short-time de- 
cay of Afoa is quartic (and not quadratic) in time, and 
is governed by the commutator (and not the difference) 
of the unperturbed and perturbed Hamiltonians, and (ii) 
for not too strong Hamiltonian perturbations, the decay 
of il/oa freezes at values inversely proportional to the 
square of the measure T of the perturbation, as defined 
by the width of the local density of states, Eq. ([TT]) . This 
allows to estimate the strength of decoherence processes 
in systems of cold trapped atoms by fitting the saturation 
value of Afoa, which is arguably easier and more precise 
than fitting decay curves over not precisely defined time 
intervals. In addition to providing an analytic deriva- 
tion of this finding, in particular relating the saturation 
level to the strength of decoherence fields, and to pre- 
dicting an initial quartic decay of AfDa, our theory gives 
an intermediate behavior of Afoa which follows that of 
the Loschmidt echo ATl- We confirmed these analytical 
findings numerically. 
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